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MASSACHUSETTS 


ABSTRACT 


A  model  is  developed  for  the  reponsivity  of  HgCdTe  detectors  employed  in  the 
GOES  Sounder  and  Imager  instruments.  Response  non-linearity  is  related  to  the 
calibrated  radiometric  accuracy  of  the  instruments.  The  model,  in  the  form  of  a  computer 
code,  could  be  utilized  to  expedite  pre-launch  calibration  of  GOES  satellite  instruments 
and  to  provide  confidence  in  GOES  in-orbit  calibration  procedures. 
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1.  INTRODUCTION 


Weather  information  to  be  provided  by  the  GOES  satellite  requires  the  satellite's 
instruments  to  be  well  calibrated.  Instrument  calibration  requires  testing  prior  to  launch  in  a 
space  simulation  chamber.  Measurements  are  taken  over  a  wide  range  of  chamber  radiometric 
conditions  and  instrument  temperatures.  Pre-launch  calibration  testing  is  relatively  time 
consuming  and  costly. 

A  model  of  the  responsivity  of  the  GOES  instruments  could  expedite  pre-launch 
calibration  testing;  however,  a  model  that  includes  instrument  response  non-linearity  has  not 
been  developed.  Non-linearity  in  GOES  instrument  response  is  attributable  to  the  HgCdTe 
detector  and  detector  preamplifier. 

This  report  describes  an  analysis  of  the  non-linear  responsivity  of  the  GOES  Sounder  and 
Imager  instruments  and  sets  forth  a  collection  of  equations,  or  model,  given  in  the  form  of  a 
computer  code  that  could  be  utilized  to  facilitate  pre-launch  calibration  of  GOES  satellite 
instruments  and  to  provide  confidence  in  GOES  in-orbit  calibration  procedures. 


1.1  BACKGROUND 


An  accuracy  goal  of  ±1K  in  scene  brightness  temperature  makes  the  GOES  instruments 
sensitive  to  small  non-linearity  in  the  relationship  between  incident  sensed  irradiance  and  the 
measured  voltage,  or  counts.  This  can  be  illustrated  by  differentiating  the  Planck  function  with 
respect  to  the  temperature,  T, 


dN  _ 
N 


dl 

T 


0) 


where  N  is  spectral  radiance,  c2=14391  fim-K,  and  X  is  the  wavelength  in  |jm.  A  plot  of  this 
expression  for  the  fractional  radiance  versus  wavelength  is  shown  in  Figure  1  for  a  ±1K 
uncertainty  in  T. 

±dN/N,% 


Figure  1.  Radiance  uncertainty  vs  wavelength  for  ±1K  uncertainty  in  300°K  scene  brightness. 
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From  Figure  1  it  is  observed  that,  particularly  at  longer  wavelengths  covered  by  the 
instruments,  the  uncertainty  in  N  must  be  small  (<  5%)  to  achieve  a  ±1K  data  accuracy.  Only 
part  of  the  radiance  uncertainty  budget  is  assigned  to  the  uncertainty  in  calibration.  Thus,  very 
small  calibration  uncertainties,  on  the  order  of  a  few  %,  are  significant. 

The  calibration  equation  converts  measured  counts  (Acounts),  taken  with  respect  to  a 
space  or  zero  radiance  reference,  to  scene  radiance.  The  equation  is  imbedded  in  NOAA 
Operational  Ground  Equipment  (OGE)  computer  software.  The  OGE  computer  software 

converts  data  to  radiance  values  from  the  Acount  data  telemetered  to  ground  from  the  orbiting 
satellite. 

In-orbit  calibration  for  GOES  is  planned  to  be  the  same  as  that  used  for  the  AVHRR: 
Each  instrument  will  be  commanded  to  sequentially  observe  space  and  an  onboard  blackbody 
target.  This  provides  data  for  determination  of  two  coefficients  in  the  calibration  equation.  For 
the  AVHRR  the  coefficients  are  the  offset  and  linear  terms. 

Consider  a  plot  of  a  calibration  equation  shown  in  Figure  2:  The  plotted  equation 
resulted  from  a  fit  to  pre-launch  calibration  data  taken  for  the  10.3  -11.3  pm  channel  of  an 
Advanced  Very  High  Resolution  Radiometer(AVHRR)  used  in  the  NOAA-7  satellitefRef  1]. 
The  instrument  employs  HgCdTe  photo-conductive  detectors  and  electronic  circuitry  that  is 
similar  to  the  circuitry  in  the  GOES  instruments.  A  non-linearity  is  evident:  The  responsivity  of 
the  AVHRR  decreases  as  the  radiance  value  increases.  To  fit  the  data  the  calibration  equation 

requires  use  of  a  constant  offset,  linear,  and  quadratic  coefficient  to  relate  radiance  to  Acounts. 


Acounts 


Figure  2.  ACounts  vs  radiance  for  channel  4  of  the  NOAA-7  AVHRR  at  a  baseplate  temperature 

of  10C. 

Coefficients  of  the  AVHRR  calibration  equation  have  been  found  to  depend  on  the 
background  flux  impinging  upon  the  detector.  Radiometric  background  flux  on  a  AVHRR 
instrument  is  influenced  by  the  instrument  baseplate  temperature.  In  Figure  2  the  baseplate 
temperature  is  10C.  At  different  baseplate  temperatures  the  calibration  curve  coefficients  are 
different  by  several  percent. 
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The  AVHRR  quadratic  coefficient  is  inferred  from  the  in-orbit  instrument  baseplate 
a  311(1  a  Pre*orblt  ground  test  data  base.  The  data  base  was  collected  by  operating  the 

AVHRR  at  various  baseplate  temperatures  in  a  space  simulation  chamber.  The  method 

SeTdi^^MtRef  If™116  AVHRR  qU3dratic  non-linearity  coefficient  is  described  by 

The  GOES  instruments  require  in-orbit  calibration  updates  over  a  range  of  fore-optics 

Certain  sJeasons  instruments  carried  aboard  the  geo-synchronous 
GOES  satellite  will  be  subject  to  radiation  heat  loading  from  the  sun.  At  such  times  the  sun 

SS  J“  Th“‘ 0f  ,,h?„fleld" of'™*  a«.d  partially  illuminates  the  instrument's  fore-optics  baffle 
change  te^omertc'S^  baffl'  ttmpera*ure  “d' thcreb^ 

Tbe  GOES  satellite  instruments  require  in-orbit  calibration  updates  over  a  range  of 

TST3?1^'  Non*hnear  instrument  responsivity  due  to  HgCdTe  detectors  will  be  a 
f  5e  detector  operating  temperature.  A  quadratic  calibration  coefficient  which 
TemStSe^  f^!.e  non-lmeanty  will  also  depend  on  the  detector  operating  temperature. 
Temperature  of  AVHRR  detectors  is  kept  constant;  however,  GOES  detectors  will  operate  at 
several  temperatures(from  94  to  105K)  over  the  course  of  a  year  and  over  the  satellite's  lifetime. 

calibration  coefficients  using  measurements  taken  in  orbit  requires  a  pre-orbit 

Er^  nf,  the  of  temPeralure  variations  expeaedin  orbU 

1  extraPolat.e  between  instrument  temperatures  could  facilitate  calibration 
equabon  coeffiaent  correction  using  a  minimum  pre-orbit  chamber  test  data  base. 


•  _  The  planned  GOES  chamber  calibration  procedure  involves  measuring  the  response  of  an 

Se MS  °f  “de?‘  n“x  ^  &<*  level  is  oblSgXSS 
the  temperature  of  a  black  body  m  the  chamber.  Allowing  the  temperature  to  stabilize  at  each 
level  consumes  considerable  chamber  test  time.  ^  31  eacn 

calihratSfem^n^rSf1  flUXJfkdS  is  dri^en  by  a  requirement  to  fit  the  measurements  to  a 
calibration  equation.  One  contribution  to  the  number  of  flux  levels  that  must  be  used  is  an 

uncertainty  in  the  form  of  the  calibration  equation.  A  model  of  the  calibration  equation  could 
remove  this  uncertainty  and,  thereby,  permit  use  of  fewer  flux  levels. 


If  Nj  is  the  number  of  flux  level  measurements,  Nc  the  number  of  coefficients,  od  the 

data  percent  precision,  and  of  the  fit  percent  uncertainty,  then  the  number  of  flux  levels  could  be 
chosen  in  accordance  with  approximately 


Vj 

crr~ - 2 — 

f  VNTTIT' 

I  C 

E.g.,  for  a  three-coefficient  calibration  equation  and  suitably  small  od,  a  minimum  of  =4-5  flux 
levels  might  be  adequate. 


12  NON-LINEARITIES 


There  are  two  causes  of  instrument  response  non-linearity: 
and  preamplifier  used  to  sense  the  changes  in  conductance.  The 


The  first  is  the  bias  circuit 
second  is  the  departure  of 
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HgCdTe  photo-conductors  from  linearity  in  response  to  impinging  radiometric  flux.  This 
distinction  is  of  interest  because  in  the  past  response  non-linearity  has  been  individually 
attributed  to  each  cause.  As  will  be  discussed,  it  is  due  to  both. 

1.2.1  CIRCUIT  NON-LINEARITY 

Figure  3  illustrates  GOES  instrument  detector  and  detector  preamplifier  circuitiy.  The 
circuit  employs  a  photo-conductive  detector  that,  along  with  a  connecting  wire  resistance, 
exhibits  a  100  to  200  Ohm  resistance,  Rd,  in  one  arm  of  a  bridge  network  The  bridge  network 
serves  as  the  input  to  a  low  noise  operational  amplifier. 

The  resistance  values  Ra  and  Rc  in  the  second  arm  of  the  bridge  are  chosen  to 
approximately  zero-balance  the  voltage  input  to  the  preamplifier  under  a  quiescent  operating 
condition. 

The  bridge  resistance  in  series  with  the  detector,  Rb,  is  used  to  control  the  detector  bias 
current  caused  by  the  bias  voltage,  Vo*  The  bias  voltage  is  —10  Volts  and  typical  bias  currents 
employed  are  *1  to  2  mA.  This  results  in  a  resistance  value  of  =5000  to  10,000  Ohms.  Because 
the  value  of  Rb  dominates  the  resistance  of  the  bridge  arm  the  influence  of  the  detector 
resistance,  Rd,  on  the  bias  current  is  small. 

The  detector's  resistance  is  a  function  of  the  power,  Pp,  due  to  incident  radiometric  flux. 
Small  changes  in  the  value  of  Rd  due  to  changes  in  Pp  cause  small  changes  in  Rd  at  the  balanced 
bridge  input  to  the  preamplifier.  These  changes  are  sensed  and  amplified  by  the  preamplifier. 


The  preamplifier  employs  a  feedback  resistor,  Rf.  The  summing  junction  in  the  feedback 
loop  is  the  negative  input  to  the  preamplifier  and  is  tied  to  node  A  shown  in  Figure  3.  The 
feedback  current.  If,  divides  at  node  A  and  returns  to  the  preamp  through  either  Ra  and  V0,  or 
through  Rc*  The  feedback  current  drives  the  voltage  at  node  A  to  be  equal  to  that  at  node  B. 
Solving  Va=Vb  gives  the  preamplifier  output  voltage,  V 
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(2) 


y  =  Y  El (^a +  Rc)  (  _  Rc  \ 

Ra^c  \Rb  +  Rd(Pp)  Ra  Rc  / 

Because  typically  Rb»Rd  and  Ra»Rc,  Eq.(2)  can  be  written  as 


Vt5(Vp)-Rd  (3) 

The  bridge  arm  resistance,  Rd,  is  the  sum  of  the  resistance  of  the  cooled  photo-conductor  and 
that  of  a  thin  connecting  wire  resistance.  The  thin  connecting  wire  in  series  with  each  contact  of 
the  detector  provides  thermal  impedance  isolation  of  the  preamplifier  circuitry  from  the 
cryogenic  cooler  that  controls  the  detector  temperature. 

If  Rw  is  the  wire  resistance  and  Rd  is  the  actual  cooled  detector  resistance,  then 

Rd  =  Rw  +  Rd  (4) 

The  value  of  Rw  is  in  the  50-75  Ohm  range  and  is  significant  with  respect  to  the  value  of  Rd. 

A  photo-conductor  exhibits  a  conductance,  Gd=l/Rd,  that  equals  the  sum  of  a  constant, 
quiescent  conductance  plus  a  photo-sensitive  conductance.  For  an  ideal  photo-conductor,  the 
photo-sensitive  conductance  is  linearly  proportional  to  impinging  radiometric  flux.  Circuit  non¬ 
linearity  in  GOES  instruments  is  attributed  to  the  fact  that  the  value  of  V,  indicated  by  Eqs.  (2) 
and  (3),  is  proportional  to  Rd  rather  than  Gd.  The  preamplifier  circuit  would  create  a  small 
response  non-linearity  even  if  an  ideal  photo-conductor  were  employed. 

1.2.2  DETECTOR  NON-LINEARITY 

Real  photo-conductors  exhibit  a  conductance^ d,  that  can  change  non-linearly  in 
proportion  to  the  incident  radiometric  flux.  Causes  of  HgCdTe  detector  photo-conductive  non¬ 
linearity  include: 

1)  non-ohmic  contacts 

2)  flux  dependence  of  the  Auger  relaxation  time* 

3)  flux  dependence  of  the  surface  recombination  relaxation  time* 

4)  flux  dependence  of  diffusion  and  electric-field  carrier  sweepout. 

The  non-linearities  considered  in  this  analysis  are  marked  with  an  asterisk. 

It  is  assumed  that  a  detector  exhibits  ohmic  contacts.  I.e.,  it  is  assumed  detector 
resistance  is  substantially  independent  of  the  direction  of  charge  flow  through  the  detector. 
During  manufacture,  detectors  are  tested  for  symmetry  of  resistance  to  forward  and  reverse  bias 
to  identify  and  eliminate  detectors  with  non-ohmic  contacts. 

Flux  dependence  of  diffusion  and  electric-field  carrier  sweepout  effects  can  be  important, 
especially  at  bias  voltages  higher  than  typically  employed  for  GOES  detectors.  However,  the 
effect  is  difficult  to  model  and  is  discussed  in  Section  2.3. 
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Flux  dependence  of  the  carrier  relaxation  times  and  carrier  sweepout  is  due  to  the  depen¬ 
dence  of  the  Auger  relaxation  time(the  principal  carrier  relaxation  mechanism),  the  diffusion 
coefficient,  and  the  ambipolar  mobility  on  carrier  density. 


6 


2.  MODEL  OF  CIRCUIT  AND  DETECTOR  NON-LINEARITY 


Detector  circuit  response  is  set  forth  as  a  series  expansion  in  a  factor  equaling  the  product 
of  the  HgCdTe  responsivity  and  the  radiometric  scene  power  incident  on  a  detector.  The 
responsivity  of  the  detector  is  then  taken  as  a  Taylor's  series  expansion  around  the  quiescent 
operating  point  determined  by  incident  background  power.  The  resulting  instrument  responsivity 
expression  is  related  to  the  calibration  equation.  Finally,  physical  causes  of  the  detector's  non¬ 
linear  responsivity  are  briefly  described. 

2.1  ANALYTIC  EXPRESSION  FOR  RESPONSIVITY 

In  orbit,  a  reference  voltage  is  determined  with  the  instrument  observing  space.  Although 
radiometric  flux  falls  on  the  detector  due  to  background  radiation  from  the  instrument,  the 
"space-look"  voltage  is  used  as  a  reference  for  zero  signal  flux  incident  on  the  instrument's 
optical  entrance  aperture.  When  the  instrument's  optical  field-of-view  is  turned  to  earth  the 
difference  voltage,AV,  between  the  "space-look"  reference  and  the  observed  voltage  is  used  to 
ascertain  the  earth  scene's  radiance. 

Data  transmitted  to  NOAA  for  processing  by  the  OGE  is  a  stream  of  digital  counts.  The 
counts  are  derived  from  an  A/D  converter  onboard  the  satellite.  It  is  assumed  that  the  A/D 
converter  is  linear  and  that  its  voltage  resolution  is  small  compared  to  the  required  precision  of 
the  data.  (The  A/D  resolution  error  contributes  a  small  random  component  to  the  instrument's 
additive  noise.) 

If  at  some  time,  t,  a  power,  Pp[t],  during  observation  of  an  earth  scene  is  incident  on  the 
detector,  then  from  Eqs.(3)  and  (4)  the  preamplifier  voltage  at  time  t  is  given  by 

V(Pp[t]  )=l^i^{RwM+Rd(Pp[t] )  -  Rett])  +  Vn[t]  (5a) 

where  Vn[t]  is  a  zero-mean  random  noise.  Similarly,  if  at  time  t+t’  a  power  Pp[t  +  1 1  from  a 
space  reference  is  incident  on  the  detector  then  the  preamplifier  voltage  at  time  t+t'  is  given  by 

V(Pp[t+f))=^-^-(Rw[t+t’]+Rd(Pp[t+t’])-  Rc[t+t’])  +  Vn[t+f]  (5b) 
The  Avoltage  reported,  AV,  is  taken  to  be 

AV  =  <  V(Pp[t])  -  V(Pp[t+f])  >  (6) 

where  the  brackets  <  >  denote  an  ensemble  average.  It  is  assumed  that  Equation  (6)  simplifies  to 
give 


AV=l^<  (R<i(pp[t) )  -  Rd(Pp[t+n ) )  > 

= 1^(<  “dCya+Pbro )  >  -  <  Rd(Pbit+ti )  >) 


(7a) 

(7b) 
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(7c) 


=  1^(Rd(Ps+Pb>'Rd(Pb)) 

where  Ps  is  the  earth  scene  power  incident  on  the  detector  and  Pb  is  the  quiescent  background 
power  from  the  instrument  environment. 

Equations  (7a)-(7c)  indicate  in  a  somewhat  formalistic  way  that  random,  temporal 
changes  due  to  noise  and  variation  in  circuit  parameters  are  taken  to  be  averaged  to  zero.  In 
practice,  data  from  both  earth  scenes  and  space  reference  scenes  can  be  averaged.  It  is  assumed 
averaging  achieves  a  reduction  of  the  random  and  temporal  noise  to  the  extent  that  it  does  not 
contribute  to  calibration  enror. 

Imprecision  of  processed  in-orbit  instrument  data  must  be  small  compared  to  the  required 
accuracy  of  the  radiometric  product.  To  be  limited  in  radiometric  accuracy  by  only  a  small 
instrument  non-linearity  is,  indeed,  a  significant  achievement. 

Equations  (7)  indicate  that  AV->  0  as  Ps  ->  0.  However,  this  has  not  been  observed  for 
the  AVHRR  or  the  SN02  and  SN03  GOES  instruments[Ref  3,4].  Chamber  test  data  indicate  a 

small  offset  I.e.,  it  was  observed  that  AV->  AV0ffset  as  Ps  ->  0  where  AVoffset  is  a  small,  but 

nonetheless  significant,  value  in  comparison  to  the  AV  associated  with  the  NEDN  and  NEDT  of 
the  instruments. 

The  AVHRR  instruments  were  tested  at  the  ITT  space  chamber  in  Ft.  Wayne,  IN.  If 
some  periodic  electromagnetic  interference  (EMI)  were  present  or  if  Pb[t]  *  Pb[t+t’]  then  such  an 
offset  might  be  anticipated.  An  offset  might  also  be  expected  if,  because  of  spurious  sources  of 
radiation  in  the  test  chamber  at  111  ,  the  reference  value  for  Ps  is  not  zero  as  it  is  assumed  to  be. 

If  an  offset  exists  due  to  the  instrumentation,  then  the  offset  should  be  included  in  the 
calibration  equation  irrespective  of  a  model.  (It  is  not  included  in  the  model  developed  from  this 
analysis.)  If,  on  the  other-hand,  the  offset  is  a  peculiarity  of  the  tests  and  is  not  exhibited  in  orbit, 
it  is  an  aberration  that  should  not  be  included  in  the  calibration  equations.  It  is  not  easy  to 
resolve  this  issue  with  the  planned  two-datum,  in-orbit  calibration  correction  procedure.  The 
offset  term  presents  an  issue  that  is  important  to  calibration;  it  is  recommended  that  it  be  further 
investigated. 

Eq.(7c)  is  in  accord  with  the  author’s  understanding  of  how  GOES  instruments  were 
designed  to  perform  in  orbit. 

The  value  of  G<j  depends  on  the  power  falling  on  the  detector  as  well  as  detector  material 
parameters: 


Gd(Pp)  =  gss  +  <lMe 


n0  + 


Vlftp)  ^(Pp)  Pp 
hcHWL 


H  W  _  1 

Le  ^d(Pp) 


(8) 


where 

Pp  =  power  impinging  on  detector.  Watts,  from  photons  at  wavelength  Xp,  pm 
gss  =  surface  state  shunt  conductance,  mho 
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q  =  charge  of  electron.  Coulomb,  h  =  Planck's  constant.  Joule-sec,  c  =  velocity  of  light, 
Mm/sec 

H  =  detector  height,  cm,  W  =  detector  width,  cm,  L  =  detector  length(along  electrodes), 
cm 

detector  photo-responsive  area  =  W  x  L,  cm2,  Le  =  detector  electrical  length,  cm. 

11  =  quantum  efficiency,  x  =  effective  bulk  relaxation  time,  sec 

(no.  Po)  =  thermally-generated  equilibrium  (electron,  hole)  concentration,  #/cm^ 

Pe  -  electron  carrier  mobility,  cm2/Volt-sec 


The  detector  parameters  gss.  M-e*  no.  t,  and  Xp  depend  on  the  detector  operating 
temperature.  For  the  sake  of  simplicity,  notation  conveying  this  dependence  will  not  be  carried. 
Substituting  Eq.(8)  into  Eq.(7c),  setting  the  gain  factor  Rf/Rc  equal  to  unity,  one  obtains  for 

Avoltage,  apart  from  an  overall  gain  factor  that  is  not  pertinent  to  the  analysis. 


AV  =  - 


P(Pb+  Ps)  PS 
!  ,  P(Pb+Ps)PS 


Rd(pb) 


where  p(Pb+ps)  is  the  detector  current  responsivity  given  by 


(9) 


I0Rd(Pb).Amp/Watt 

and  the  bias  cunent  Iq  is  given  by  V<>/Rb. 


(10) 


Equation  (9)  indicates  that  AV  is  proportional  to  the  scene  power  incident  on  the  detector, 
Ps,  apart  from  a  term  in  the  denominator  of  the  right-hand-side  that  is  also  proportional  to  Ps. 
This  denominator  term  produces  a  small  non-linearity  in  the  response  that  is  attributable  to  the 
circuit  non-linearity. 

For  small  non-linearity  Eq.(9)  can  be  approximated  by 


AV  =  -p(Pb+Ps)  Rd(Pb)Ps  1- 


Pgb+ps)ps^PS)j2ps2j, 


(11) 


Equation  (11)  indicates  that  the  circuit  non-linearity  may  be  accounted  for  by  a  polynomial 
correction  factor  in  P §.  It  is  important  to  note  that  this  does  not  include  a  detector  response  non¬ 
linearity.  The  detector  non-linearity  lies  in  the  detector  current  responsivity,  p(Pb+Ps).  If  the 
detector  response  non-linearity  is  small  it  can  be  approximated  by  the  first  terms  of  a  Taylor's 
Series 


p(Pb+Ps)  =  p(Pb)  +  p'(Pb)Ps  +  0.5  p"  (Pt>)Ps2  + .  (12) 

which  is  a  reasonable  approximation  for  Ps  «  Pb.  Substituting  Eq.(12)  into  Eq.(9)  gives  an 
expression  that  includes  both  the  circuit  and  detector  responsivity  non-linearities: 
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AV  =  -p(Pb)  Rd(Pb)Ps 


(13) 


1- 

IpCPb) 

p’(Pb)l 

P  _ 

fp(Pb)\2  2p'(Pb)  p"(Pb)| 

p  2 

l  Io 

p(Pb)l 

rs 

l  Io  J  Io  P(Pb)  I 

rs 

Equation  (13)  is  a  simple  expression  for  a  response  non-linearity  in  Ps  up  to  cubic  order. 
Note  there  are  two  contributions  to  the  coefficients  of  the  non-linear  terms:  One  is  traceable  to 
the  circuit  non-linearity  and  one  is  attributable  to  the  detector  response  non-linearity.  In  the 

cubic  coefficient  the  effects  are  cross-coupled.  Since  p'(Pb)  is  negative  (p  decreases  with 
increasing  Pb)  the  circuit  and  detector  non-linearity  effectively  add  together  in  the  coefficient  for 
the  quadratic  term. 

2.2  INSTRUMENT  RESPONSIVTTY  AND  THE  CALIBRATION  EQUATION 

It  is  useful  to  examine  how  the  responsivity  of  the  instrument,  Eq.(13),  relates  to  a  non¬ 
linear  calibration  equation  used  by  NOAA  in  the  OGE. 

Apart  from  an  overall  electronics  gain  factor,  which  is  set  to  unity  for  this  analysis,  the 
AV  responsivity  of  the  instrument  to  power  on  the  detector,  Ps,  can  be  expressed  as 

AV  =  ocj  P$  +  0-2  Ps2  +  ct3  Ps3  + .  (14a) 


where  the  coefficients  (aj,  aj, . )  are  taken  from  Eq.(13).  A  calibration  equation  in  terms  of 

AV  and  power  on  the  detector  can  be  expressed  in  the  form 

Ps  =  Pi  AV  +  p2Av2  +  P3Ay3  + .  (14b) 

Substituting  Eq.(14a)  into  (14b)  and  solving  for  the  relation  between  the  coefficients 
gives,  to  third  order  in  AV  and  Ps 


<*2 

(ai)3 


and  P3 = 


2(02)2 

-  “1 


-  <x3 


(01)4 


(15) 


Equation  (15)  indicates  that,  for  small  corrections  about  a  linear  relationship,  one  may  examine 
the  non-linearity  of  calibration  either  in  the  instrument  responsivity  formulation  or  the  calibration 
equation  formulation.  Results  may  be  mapped  from  one  to  the  other  via  relationships  of  the 
coefficients  of  a  series  expansion  in  either  formulation. 


The  OGE  converts  Acounts  reported  from  an  instrument  channel  into  an  observed  spectral 
radiance.  To  describe  this  process  consider  two  variables  (x,y)  defined  as 


Pss  10 


,-3 


(VVoAv)y=V’  y 


in 


mW 


nr  -  Sr-cm 


(16a) 


where 

y  =  instrument  channel  measured  spectral  radiance 
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Kp  =  instrument  channel  spectral  radiance  to  detector  scene  power  conversion  factor 

Ao  =  instrument  entrance  aperture  area,  m2 

Qo  =  instrument  IFOV,  Sr 

Av  =  instrument  effective  bandpass,  cm'l 

to  =  instrument  optical  transmission 

and 

AV  s  (a  Va^)  x  =  kv  x  (  xinAcounts  (16b) 

where 

x  =  instrument  channel  reported  Acounts 

Kv  =  AVa/d,  instrument  A/D  (Avoltage-to-counts)  conversion  ratio 


The  calibration  equation  used  by  the  OGE  is  of  the  form 

y  =  bi  x  +  l>2  x2  +  b3  x3+ .  (17) 

Utilizing  Eqs.(16a)  and  (16b)  and  comparing  Eq.(17)  with  Eq.(14b)  it  is  found  that 


bl=^Pl>  b2=  P2  and  b3=  ^-03  (18) 

Rp  Kp  Kp 

It  is  possible  to  simply  relate  the  OGE  calibration  equation  coefficients  to  the  coefficients  of 
Eq.(14b).  It  is  notationally  convenient  to  minimize  extraneous  parameters  that  vary  from 

channel-to-channel.  (The  instrument  electronic  gain  is  incorporated  in  the  coefficient  kv.)  In 
this  analysis  the  P  coefficients  of  Eq.  (14b)  are  examined. 


From  Eqs.  (13),  (16)  and  (18),  the  coefficients  bi  and  b2  are  given  by 


bl  =  - 


4VVd  1 

1 

.lO-3  (ao£20Co  Av)_ 

PRd 

b2=  + 


. . ""I 

4 

_ 1 

[p(Pb)  p'(Pb)l 

1 

.  10-3  (AqQqIq  Av)J 

[  Io  P(Pb)J 

(pRd)2 

(19a) 


(19b) 


Equations  (19a)  and  (19b)  suggest  that  the  dependence  of  the  coefficients  in  the  OGE 
calibration  equation  on  background  power  and  temperature  can  be  assessed  from  detector  and 

detector-preamplifier  testing.  For  example,  the  product  (p  Rd)  is  the  voltage  responsivity  in 
Volts/Watt.  This  quantity  is  measured  to  ascertain  the  suitability  of  manufactured  detectors  for 
meeting  specification.  The  detector  resistance  is  also  measured  using  low  resistance  connecting 
wires.  These  measurements  are  routinely  carried  out  under  controlled  detector  current  bias,  I0,  a 
single  background  power,  Pb,  and  at  several  detector  temperature  conditions.  The  change  of 
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detector  responsivity  with  background,  p'(Pb),  can  be  determined  by  measuring  the  responsivity 
at  two  or  more  controlled  background  flux  levels. 

Present  measurements  might  be  extended  over  anticipated  instrument  operating 
conditions  in  bias  current,  detector  temperature,  and  background  power.  The  results  could  be 
analyzed  using  Eqs.  (19)  to  indicate  the  effects  of  the  detector  and  preamplifier  circuit 
performance  on  the  calibration  equation  before  selection  of  detectors  and  assembly  into  an 
instrument.  (The  test  dewar  and  set-up  would  have  to  accurately  reproduce  the  background  flux, 
signal  flux,  and  temperature  that  the  detectors  experience  in-orbit.) 


2.3  DETECTOR  RESPONSE  NON-LINEARITY 

The  detector  response  non-linearity  is  primarily  due  to  the  influence  of  a  photo-excited 

carrier  density  on  the  re-combination  or  relaxation  time,! ,  of  photo-excited  carriers.  Equations 
that  describe  detector  physics  are  given  in  Appendix  A. 

The  major  factor  contributing  to  HgCdTe  response  non-linearity  is  dependence  of  the 
Auger  relaxation  time.XA, on  carrier  density 


x  =  !AAi  = 


2(ni)^Ai 


(no  +  p0+np)(no  +  np) 


Al 


(20) 


where 

nj  =  intrinsic  electron  density,  cm' 3 
no  =  equilibrium  electron  density,  cm' 3 
p0  =  equilibrium  hole  density,  cm'3 
np  =  photo-excited  carrier  density,  cm'3 

TAi =  intrinsic  Auger  relaxation  time(np  =  0  and  no=po=ni)>  Sec 
Ai  =  Aj(np)  surface  recombination  reduction  factor. 

For  a  steady-state  incident  radiometric  power,  Pp,  the  steady-state  photo-excited  carrier  density, 
np,  is  given  by 


TJ  Xp  T  Pp 

"P - teHAj 


(21) 


The  detector's  photo-sensitive  conductivity  is  proportional  to  np  which,  through  Eq.(21), 
is  also  proportional  to  Pp.  However,  from  Eq.(20)  it  is  seen  that,  as  np  approaches  a  density 
nearing  no,  the  relaxation  time  will  tend  to  decrease.  Through  Eq.(21)  the  decreasing  relaxation 
time  decreases  the  responsivity  of  the  photo-sensitive  conductivity  to  Pp.  This  leads  to  a  non¬ 
linearity  in  the  response  to  Pp. 


The  factor  Ai  is  dependent  on  np  through  the  carrier-dependent  diffusion  coefficient  of 
HgCdTe  material  and  also  influences  the  non-linearity,  although  less  significantly  than  x^. 
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Diffusion  and  electric-field  induced  carrier  sweepout  also  affect  the  responsivity  of  the 

photo- sensitive  conductivity  to  Pp.  These  mechanisms  influence  responsivity,  p(Ps),  as  a 
function  of  the  position  between  detector  electrodes  and  tend  to  reduce  the  small-signal,  linear 
responsivity  of  a  detector.  The  reduction  of  the  small-signal,  linear  responsivity  due  to  these 
mechanisms  is  described  in  an  appendix  of  the  NOAA  Detector  Team  Final  Report[Ref  5].  To 
the  author's  knowledge,  the  influence  of  these  mechanisms  on  the  non-linearity  of  the  detector 
response  has  not  been  reported  in  the  literature. 

Modeling  the  responsivity  of  a  detector  affected  by  diffusion  and  electric  field  sweepout 
requires  considering  np->np(x),  tA(np)*>xA(np(x))*  and  p(Ps)->p(x,Ps),where  x  is  a  distance 
along  the  axis  between  electrodes.  This  is  difficult  to  analyze  and  model;  however, 
experimentally  measuring  p(x,Ps)  provides  an  evaluation  of  diffusion  and  electric-field  induced 
carrier  sweepout  effects.  A  uniform  p(x,Ps)  indicates  absence  of  effects.  LIRIS  measurements 

indicate  pfxj’s)  becomes  just  noticeably  non-uniform  at  bias  currents  of  2.0-3.0  mA  or  larger. 
At  a  bias  current  of  1.0  mA,  no  spatial  non-uniformity  could  be  discemed[  Ref  6]. 

The  analysis  ignores  diffusion  and  electric-field  induced  carrier  sweepout.  It  should  be 
applicable  to  GOES  detectors  having  bias  currents  up  to  =2.0-  3.0  mA.  Pursuit  of  the  effects  of 
diffusion  and  electric-field  induced  carrier  sweepout  on  response  non-linearity  is  recommended 
as  a  useful  addition  to  the  model. 
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3.  DESCRIPTION  OF  MODEL  CALCULATIONS 


A  computational  method  used  to  obtain  a  solution  for  the  response  is  described.  The 
solution  gives  a  responsivity,  9l(Pb,Ps).  as  a  function  of  Ps.  The  percent  change  in  the  small- 

signal  responsivity,  9i(Pb,0),  and  9UPbJPs)  is  defined  as  a  percent  non-linearity  function 
%NL(Pb,Ps)-  The  computed  quantities  are  related  to  the  calibration  coefficient  parameters 
described  in  Section  2.  A  brief  review  is  given  of  detector  data  available  to  the  author  at  the  time 
the  analysis  was  undertaken. 

3.1  ITERATIVE  SOLUTION 

Considering  Ala  constant  and  independently  determining  no,  po,  ni,  TAi,  Ti,  and  Xp  for  a 
detector  allows  one  to  substitute  Eq.(20)  into  Eq.(21)  and  obtain  a  cubic  equation  for  np  in  terms 
of  Pp.  This  equation  has  only  one  real  root  and  produces  a  solution  that  is  linear  in  Pp  for 
np/no«l.  The  task  of  the  iteration  is  to  find  np  that  yields  a  self-consistent  solution  to  Eqs.(20) 
and  (21). 

The  iteration  commences  by  first  setting  Ai  =  1,  the  value  of  Ai  in  the  absence  of  surface 
recombination  effects.  Then,  a  zeroth-order  value  of  np  is  determined  from  Eq.(21)  using  a  value 

of  x  derived  from  Eq.(20)  with  np  =  0.  This  zeroth-order  value  of  np,  np[0],  is  used  to  calculate 
Ai(np[0]).  The  value  of  Ai(np[t)])  is  then  treated  as  a  constant  and  used  in  the  cubic  equation 
that  results  from  the  simultaneous  solution  of  Eqs.(20)  and  (21)  just  described.  The  resultant 

value  of  np,  np[l],  is  used  to  calculate  t(np[l])  and  Ai(np[l]).  These,  in  turn,  are  used  to  modify 
the  cubic  equation  which  is  then  solved  for  an  improved  estimate,  np[2]. 


This  process  is  repeated  until  the  difference  between  the  estimate  of  np[k]  and  np[k+l] 
becomes  suitably  small. 


3.2  CALIBRATION  EQUATION  COEFFICIENTS 


The  iterative  procedure  is  used  to  compute  a  voltage,  Vb,  that  is  proportional  to  the 
preamplifier  voltage  produced  by  the  background  power  in  accordance  with  Eq.  (2).  This 
calculation  ignores  the  gain  factor  and  bridge  offset  voltage  subtracted  from  Vb  given  in  Eq.(2). 


Rd(Ph)  \ 

Rb  +  Rd(Pb)  / 


(22a) 


A  voltage  given  by  the  background  power  plus  signal  power,  Vbs,  is  computed  in  a  similar 
manner. 


Vbs  =  V0 


/  RdOVPs)  -l 
lRb  +  Rd(Pb+Ps)/ 


(22b) 
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A  difference  voltage,  AVbs,  is  evaluated  that  gives  the  difference  between  Vbs  and  Vb  and 
represents  the  AV  reported  by  the  instrument  in  accordance  with  Eq.(7c). 


AVbs  =  Vbs  *Vb 


The  voltage,  AVbs.  gives  a  solution  from  the  iterative  procedure  for  the  response  of  an 
instrument  to  an  incident  scene  power,  Ps,  for  a  given  Pb-  The  voltage  responsivity  of  the 

instrument  is  taken  to  be  91  (Pb,Ps)  where 


03) 

*s 

The  response  non-linearity  for  a  given  Pb  can  be  envisioned  as  the  deviation  of  91  (PbJ*s)  from  a 
straight  line  when  plotted  as  a  function  of  Ps.  It  is  difficult  to  quantitatively  determine  small 
non-linearities  from  such  a  plot.  For  this  reason  the  computations  include  a  percent  non-linearity 

function,  %NL(PbJPs).  that  compares  91  (Pb J*s)  with  91  (Pb.0). 

To  calculate  9t(Pb,0)  a  small-signal  voltage,  Vbss.  is  calculated  as  in  Eq.  (22b)  where 
Ps->Pss=10-10Pb. 


Vbss  ~  V0 


A  AVbss  is  then  computed 


AVbss  =  Vbss -Vb 


and  the  small-signal  responsivity,  91  (Pb,0),  is  evaluated  as 

9l(Pb»0)  =  ^kss..  (24) 

pss 

The  deviation  of  9t(Pb4>s)  from  9l(Pb,0)  is  a  measure  of  the  departure  of  the  responsivity 
from  linearity.  The  percent  non-linearity  function  is  defined  as 

%NL(Pb.Ps>  =  100  foPW  .  l)  (25) 

\9t(Pb,0)  I 

The  percent  non-linearity  can  be  interpreted  in  the  following  way:  It  is  the  percent  reduction  in 
actual  AVbs  for  Pb+Ps  compared  to  the  AVbs  that  would  occur  for  Ps  due  to  the  small-signal 
responsivity  (Ps«Pb)  value. 
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Comparing  the  results  of  the  calculations  with  the  analytic  expression  of  Eq.(13),  gives 
the  correspondences: 


SR(Pb,0)  =  -p(Pb)Rd(Pb) 


(26) 


and 


%NL(Pb,Ps)  s  -100 


p(Pb)  P'Wlp  |/p(Pb)f ,  2p'(pb)  05P,,(pb)lp  2 
17"pw)s'lio)  Io  P(Pb)|SJ 


(27) 


Equation  (27)  relates  the  %NL  to  p(Pp)  and  detector  current,  Io- 


The  percent  non-linearity  function  and  SR  (Pb,0)  are  used  to  evaluate  the  effect  of  the 
responsivity  on  the  coefficients  of  the  calibration  equation.  Comparing  the  calculated  quantities 
to  Eq.(14a)  leads  to  the  identification 

ai=9t(Pb.0>  and  =  (28) 


Using  Equation  (15),  the  values  of  pj  and  P2  are  evaluated  to  examine  the  change  of  the 
calibration  equation  coefficients  with  temperature  and  background  power. 

The  non-linearity  evaluated  by  the  preceding  computations  includes  both  circuit  and 
detector  non-linearities.  The  detector  non-linearity,  however,  is  due  solely  to  the  change  in  the 

carrier  relaxation  time,  t.  To  distinguish  between  circuit  and  detector  non-linearity  a  %NLt 
function  is  evaluated: 


/t(P.+Pc)  \ 

(29) 

The  %NLx  function  gives  the  contribution  of  the  detector  responsivity  non-linearity  to  the 
overall  %NL. 

3.3  AVAILABLE  DETECTOR  DATA 

The  HgCdTe  detectors  used  in  GOES  are  tested  and  characterized  by  the  manufacturer 
and,  independently,  by  ITT  before  being  mounted  in  a  an  instrument  package.  The  tests  are 
currently  conducted  at  a  single  background  flux  level  and  do  not  characterize  the  detector's  non¬ 
linearity. 

The  values  of  no,  Po.  ni,  TAi,  T|,  and  Xp  depend  sharply  on  the  temperature,  extrinsic 
donor  impurity  concentration,  N<j,  and  the  molar  fraction,  x,  in  the  Hg(l-x)Cd(x)Te  material.  To 
assess  approximate  values  for  these  parameters  a  match  was  empirically  made  to  data  collected 
on  the  detectors  by  Ed  Andrews  of  ITT[Ref  7].  Additional  data  have  been  supplied  by  Jean 
Hartley  [Ref  8], 
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4.  EXAMPLES  OF  MODEL  CALCULATIONS 

Two  examples  of  model  calculations  are  given  for  detectors  to  be  used  in  the  GOES 
SN03  satellite.  Estimating  detector  parameters  from  test  data  is  the  first  step.  Variation  of  the 
small-signal  responsivity  and  percent  non-linearity,  %NL,  versus  background  power  and  detector 

temperature  is  evaluated.  Finally,  variation  of  the  calibration  equation  coefficients  Pi  and  p2 
versus  background  power  and  detector  temperature  is  calculated. 


4.1  FITTING  PARAMETERS 

Parameters  of  the  model  are  derived  from  detector  test  data.  As  described  in  Appendix 
A,  one  of  the  more  critical  parameters  is  the  Cd  molar  fraction,  x,  in  the  Hg(l-x)Cd(x)Te  alloy. 
This  can  be  deduced  from  spectral  measurements  of  the  relative  responsivity.  This  is  illustrated 
in  Figure  4  and  Figure  5. 

The  50%  relative  responsivity  wavelength  in  the  cut-off  region  is  a  sensitive  function  of 
the  detector  operating  temperature.  By  comparing  the  observed  50%  relative  responsivity  point, 
shown  as  a  black  dot,  to  the  relative  responsivity  predicted  by  the  equations  for  the  responsivity 
given  in  Appendix  A,  a  value  of  x  can  be  determined  that  matches  the  data. 


Relative  Responsivity 


Wvl,  fim 


12  14  16  18 

Figure  4.  Match  of  50%  relative  responsivity  wavelength  vs  temperature  for  the  Judson 

LongWave  J-033  detector. 


For  the  case  of  the  LIRIS  I3-L009  detector  the  spectral  response  of  the  LIRIS  I3-L006 
has  been  used.  Spectral  response  data  on  L009  was  not  available,  but  was  stated  to  be  virtually 
identical  to  L006  [Ref  9].  } 
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Relative  Responsivity 
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Figure  5.  Match  of  50%  relative  responsivity  wavelength  vs  temperature  for  the  LIRIS I3-L006 

detector. 


It  is  seen  from  Figures  4  and  5  that  the  values  of  x  are  0.1956  and  0.2338,  respectively. 
These  two  values  are  for  the  smallest  and  largest  bandgap  material  employed  in  the  SN03 
instrument,  respectively. 

The  model  equations  given  in  Appendix  A  have  a  precision  that  is  sensitive  to  changes  in 
the  fourth  decimal  place  of  x.  However,  the  absolute  accuracy  of  the  equations  and  the  predicted 
value  of  x  is  on  the  order  of  =  ±2%. 

Matches  to  other  parameters  were  obtained  by  comparing  the  model  to  measured 
responsivity  and  detector  resistance  versus  bias  current.  The  results  of  this  match  is  given  in 
Table  1.  The  optical  area  and  bias  current  assumed  is  also  included  in  Table  1. 

TABLE  1 


Detector  Parameters 


Pararneigr 

I3-LQQ2 

J-033 

Cd  molar  fraction,  x 

0.2338 

0.1956 

optical  area,  cm^ 

1.07  x  IO"4 

3.77  x  IO’5 

Nd,  #/cm3 

5  x  1014 

3  x  1014 

S,  cm/sec 

50 

50 

gss,  rnMho 

4 

7 

Le/L 

4 

3 

Rd,  Ohms 

=200(1 05K) 

=90(1 02K) 

SR[Pb,  0],  kV/Watt 

=120(105K) 

42(1 02K) 

Io,  mA 

1 

1 
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4.2.  SMALL-SIGNAL  RESPONSIVITY  VERSUS  BACKGROUND  AND 
TEMPERATURE 

Variation  of  instrument  small-signal  responsivity  with  changes  in  background  power  and 
detector  temperature  is  a  critical  issuefRef  10].  To  evaluate  the  responsivity  from  the  models  it 
is  necessary  to  estimate  the  background  power  and  maximum  scene  power  falling  on  the 
detectors.  Table  2  is  a  summary  of  channel  characteristics  assumed. 

TABLE  2 


Channel  Parameters 


» 

Parameter 

Imager  3 

LW  Sounder  Ch  7 

<X> 

6.75  pm 

12  pm 

$ 

AX 

0.5  pm 

0.7  pm 

Max  TScene 

320K 

330K 

Pb 

12  nWatts 

110  nWatts 

Max  Ps 

7  nWatts 

16  nWatts 

The  results  of  model  computations  are  shown  in  Figures  6  and  7.  Note  that  in  the  Figures 
the  responsivity  in  kV/  Watt  is  a  negative  number. 


It  is  seen  from  Figure  6  that  both  the  Imager  Channel  3  and  the  Sounder  Channel  7  exhibit 
a  small-signal  responsivity  that  decreases  with  increasing  background  power.  The  decrease  is 

attributable  to  the  decrease  in  earner  relaxation  time,t,  caused  by  the  increasing  photo-excited 
carrier  density,  np,  that  results  from  increasing  Pt>.  This  effect  is  indicated  by  Eqs.(20)  and  (21). 
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Figure  7.  Variation  of  the  small  -signal  responsivity  of  Imager  3  and  Long  Wave  Sounder 
Channel  7  vs  detector  temperature.  ( The  background  power  on  the  detector  is  21  nW for  the 

Imager  and  110  nW for  the  Sounder.) 

Figure  7  shows  a  difference  between  the  small-signal  responsivity  versus  detector 
temperature  for  the  small  and  large  bandgap  materials.  From  Eq.(8)  it  can  be  reckoned  that  the 
temperature  dependence  of  the  small-signal  responsivity  is  influenced  by  three  effects:  The 

electron  mobility,  |ie,  the  carrier  relaxation  time,  x,  and  the  equilibrium  electron  density,  no, 
which  tends  to  shunt  changes  in  photo-excited  conductance. 

Electron  mobility  increases  with  decreasing  temperature,  tending  to  increase  the 
responsivity  of  both  detectors  with  decreasing  temperature.  However,  the  carrier  relaxation  time 
and  equilibrium  electron  density  dominate  the  temperature  dependence  over  the  operating  range 
shown  in  Figure  7.  These  two  dominant  effects  behave  differently  for  the  two  detectors  in  the 
operating  temperature  range.  To  heuristically  explain  this  difference  in  behavior,  reference  is 
made  to  equations  in  Appendix  A. 

The  Sounder  Channel  7  detector  has  a  valence-to-conduction  bandgap  that  is  relatively 
small  compared  with  its  thermal  energy  over  the  operating  temperature  range.  The  small 
bandgap  causes  the  value  of  nj  to  be  large  and,  in  particular,  large  compared  to  Nd-  This,  in  turn, 
causes  n©  and  po  to  equal,  or  nearly  equal,  the  intrinsic  carrier  density,  nj.  The  intrinsic  carrier 
density,  ni,  decreases  exponentially  with  decreasing  temperature.  This  carrier  concentration 
effect  is  indicated  in  Eqs.  (A5)  and  (A6). 

From  Eq.(20)  or  Eq.(A9b)  it  can  be  seen  that  for  no=Po~ni  the  carrier  density  factors  in 
the  numerator  and  denominator  tend  to  cancel,  making  the  carrier  relaxation  time  proportional  to 

the  intrinsic  Auger  relaxation  time,  T^i-  The  intrinsic  Auger  relaxation  time,  T^j,  increases  with 
decreasing  temperature  as  indicated  in  Eq.(A9a). 

The  increase  of  the  responsivity  of  the  Sounder  Channel  7  detector  with  decreasing 
temperature  results  from  the  decrease  of  no(less  shunt  conductance)  and  the  increase  of 

x^i(longer  carrier  lifetime)  with  decreasing  temperature. 

The  Imager  3  detector  has  a  bandgap  that  is  relatively  large  compared  with  its  thermal 
energy  over  the  operating  temperature  range.  This  causes  the  value  of  nj  to  be  relatively  small 
and  on  the  order  of  Nd  over  the  operating  temperature  range.  This  in  turn,  causes  n0  to  equal,  or 
nearly  equal,  the  detector’s  donor  impurity  density,  Nd,  and  the  value  of  p0  to  be  small.  This 
carrier  concentration  effect  also  results  from  Eqs.  (A5)  and  (A6).  From  Eq.(20)  or  Eq.(A9b)  it 
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can  be  seen  that  for  Po«n0=N(i  the  carrier  relaxation  time  is  proportional  to  the  product,  ni2!^. 

The  product  ni2x^j,  degrease?  with  decreasing  temperature  due  to  the  strong,  exponential 
dependence  of  nf  on  temperature  as  given  by  the  bj  factor  in  Eq.(A5). 

hlgheT  ?Perating  temperatures,  decreasing  the  temperature  increases  the  Imager  3 

iS?LnS?!nSIVlty'  — !S  1S  btcaus.e  *e  mobility  increases  and  because  the  detector  Kin  a 
earner  density  regime  similar  to  that  of  the  Sounder  Channel  7  detector  described  above. 

Upon  decreasing  temperature,  a  temperature  is  reached  for  which  no=N<j.  At  this 
temperature  the  earner  relaxation  time,  x,  becomes  proportional  to  the  product  ni2  XAi  The 

SfiSj deT“  *?‘h  dK?eaf,8  temperature  due  to  the  decrease  of 
SSed  in  h|uiI  7  '  Weak  maxunum  In  the  Imag«  3  small-signal  responsivity 

•  •  .J^e  model  demonstrates  that  the  temperature  dependence  of  the  small-signal  responsivitv 

2* hSvS  oS  m  KU“Ced  by  J1*  d°Ping  level>  Nd,  the  Cd  molar  fraction,  x,  vSich  determine 
the  HgCdTe  alloy  bandgap,  and  by  the  operating  temperature  range. 


43  PERCENT  NON-LINEARITY 

'sttttwte  of  the  percent  non-linearity  is  calculated  for  the  two  detectors.  Results  of  the 
calculations  are  shown  in  Figure  8.  Note  that  the  value  of  %NL  shown  is  negative  TTie  Dement 
non-lmeanty  for  a  fixed  background  flux  increases  with  increasing  signal  or  feene  power 

.  ,  Th®  Percent  non-linearity  reaches  =4-6%  at  high  scene  flux  levels.  Based  on  the 
inference  drawn  from  Eq.(l)  and  from  Figure  1,  this  is  a  significant  non-linearitv  at  loneer 

k  ane!mi?SS  T*  regard  to  cajjtffation.  Thus,  the  detector  and  preamplifier  circuit  non-lineanty 
is  an  important  issue  to  the  GOES  instrument  calibration.  y 

%NL  vSiKP^wi^'SnWfrh  Eq-(15  indiucatf  that  a  qu^atic  or  higher-order  variation  in 
ff  ?s  u  duiilTfrom  a  cublc  or  hlgher-order  dependence  of  the  responsivity  on  Pc 
Figure  8  shows  that  the  %NL  is  linear  with  respect  to  Ps,  or  nearly  so,  and  dJes  not  exhTbit  a 

significant  higher-order  dependence  on  Ps.  Thus,  the  responsivity,  9t(Pb,Ps),  is  suitably 

represented  by  a  quadratic  function  in  Ps  for  a  given  background  flux  level.  °  y 


%NL 


Imager  3 


%NL 


Sounder  7 


nWatt 


-nWatt 


H  °  0  ™  *  10  15  20  "25 

Flhnnt 8'  Per,cent  nonjlnearity  vs  signal  power,  Ps.  (Operating  temperature  is  100K  with 
background  power  ofldnW  and  llOnW  on  13  and  the  LW  Sounder  Ch  7,  respectively  ) 
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The  percent  non-linearity  associated  with  the  detector,  %NLx,  given  in  Eq.(29)  was 
evaluated.  It  was  found  that  approximately  70%  of  the  non-linearity  was  due  to  the  dependence 
of  the  carrier  relaxation  time  on  the  photo-excited  earner  density,  np,  as  given  by  Eq.(20).  I.e., 
for  the  SN03  Imager  3  and  Sounder  Channel  7  the  detector  non-linearity,  rather  than  the  circuit 
non-linearity,  was  found  to  cause  =70%  of  the  overall  non-linearity. 


4.4  CALIBRATION  EQUATION  COEFFICIENTS 

The  model  calibration  equation  is  a  polynomial  in  Acounts  as  desenbed  in  Section  2.2. 
For  a  quadratic  fit,  the  coefficients  of  the  Acount  terms,  bi  and  b2,  are  proportional  to  parameters 

pi  and  p2,  respectively,  as  given  by  Eq(18). 

In-orbit  correction  of  calibration  coefficients  relies  on  two  data  points  collected  for  each 
channel.  The  data  points  are  to  be  used  to  adjust  offset  term(not  in  the  model)  and  the  linear 

term.  Pi.  The  quadratic  term,  P2,  is  to  be  adjusted  using  telemetered  temperature  data  from  the 
satellite.  Predicted  in-orbit  variation  of  pi  and  p2  due  to  detector  and  detector  circuit 
responsivity  is  a  goal  of  the  modeling  effort. 


Figures  9  and  10  show  the  variation  of  Pi  versus  background  power  and  detector 
temperature.  Note  the  values  of  P  are  negative. 


Betal*l(T6 
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Sounder  7 
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Figure  9.  Variation  of  pi  for  Imager  3  and  Long  Wave  Sounder  Channel  7  vs  background 
power.  (The  temperature  of  the  detectors  is  100K.) 


Figure  9  shows  that  for  both  Imager  3  and  Long  Wave  Sounder  Channel  7  pi  increases  in 
magnitude  with  increasing  background  power.  This  is  consistent  with  the  dependence  of  the 
small-signal  responsivity  versus  background  power  shown  in  Figure  6:  The  responsivity 
decreases  in  magnitude,  causing  the  calibrated  estimate  of  impinging  radiance  given  by  the 
calibration  equation  to  increase  for  a  given  value  of  Acounts. 
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Figure  10.  Variation  of  f}]  of  Imager  3  and  Long  Wave  Sounder  Channel  7  vs  detector 
temperature.  (The  background  power  on  the  detector  is  21  nW  for  the  Imager  and  110  nWfor 

the  Sounder.) 


The  variation  of  pi  to  changes  in  detector  temperature  for  Imager  3  and  Sounder  Channel 
7  shown  in  Figure  10  reflect  the  comparative  small-signal  reponsivity  versus  detector 
temperature  shown  in  Figure  7:  Imager  3  exhibits  a  local  minimum  in  pi.  For  Sounder  Channel 
7  the  magnitude  of  Pi  increases  with  increasing  temperature  and  is  more  sensitive  to  temperature 
than  the  Pi  coefficient  of  Imager  3. 


Figures  11  and  12  show  the  effects  of  background  power  and  temperature  predicted  by 

the  model  on  the  p2  coefficient.  This  coefficient  is  to  be  corrected  using  in-orbit  calibration 
measurements  and  chamber  test  data  as  described  in  Section  1.1. 


Beta2*10~23  Beta2*l(T21 

Imager  3  Sounder  7 


Figure  11.  Variation  of  fy  of  Imager  3  and  Long  Wave  Sounder  Channel  7  versus  background 
power.  (The  temperature  of  the  detectors  is  100K.) 


Figure  11  indicates  that  the  p2  coefficient  increases  with  increasing  background  power 

for  both  Imager  3  and  Sounder  Channel  7.  The  p2  coefficient  of  the  Sounder  Channel  7  is 
significantly  more  sensitive  to  background  power  than  the  Imager  3  channel. 
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Figure  12.  Variation  of  @2  of  Imager  3  and  Long  Wave  Sounder  Channel  7  vs  detector 
temperature.  ( The  background  power  on  the  detector  is  21  nW  for  the  Imager  and  110  nW  for 

the  Sounder.) 

Figure  12  indicates  that  the  variation  of  P2  with  detector  temperature  decreases  with 

increasing  temperature  over  the  operating  range  while  the  P2  coefficient  of  the  Sounder  Channel 
7  increases.  The  sensitivity  of  the  Sounder  Channel  7  coefficient  to  temperature  is  significantly 
greater  than  Imager  3. 

Figures  11  and  12  indicate  that  the  variation  of  02  with  background  power(optics  baffle 
temperature)  and  detector  temperature  depends  in  an  approximately  quadratic  manner  on  these 
parameters.  Interpolation  between  pre-orbit  chamber  test  data  points  and  in-orbit  measurements 

will  be  required  to  correct  f)2-  The  model  indicates  that  a  quadratic  interpolation  scheme  should 
be  satisfactory. 
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5.  SUMMARY 


A  small,  few  percent(<  5%)  non-linearity  in  responsivity,  if  not  acounted  for  in 
instrument  calibration,  could  significantly  affect  the  accuracy  of  the  products  provided  by  the 
GOES  satellite,  especially  at  longer  wavelengths.  Determining  response  non-linearity  is  a  driver 
in  the  amount  of  required  pre-launch  instrument  chamber  testing.  Response  non-linearity  is  due 
to  the  HgCdTe  detector  and  preamplifier  circuitry  employed  by  GOES  instruments. 

Analysis  of  HgCdTe  detector  and  preamplifier  response  non-linearity  was  undertaken 
based  on  detector  physics  reported  in  the  literature.  Although  the  analysis  did  not  consider  flux- 
dependent  non-linearity  resulting  from  sweepout  effects  due  diffusion  and  electric  bias  fields,  it 
should  apply  to  GOES  instruments  employing  detectors  biased  up  to  =2.0-3.0  mA. 

The  analysis  also  related  instrument  responsivity  and  the  OGE  calibration  equation.  No 
rationale  was  uncovered  for  a  constant  offset  term  in  the  OGE  calibration  equation,  although 
such  a  term  is  employed  in  practice  based  on  pre-launch  chamber  tests.  Further  pursuit  of  the 
cause  of  the  offset  term  is  recommended. 

A  computer  code  model  of  the  GOES  instrument  responsivity  was  developed.  The 
objective  was  to  develop  a  model  that  could  be  utilized  to  expedite  pre-launch  calibration  of 
GOES  satellite  instruments  and  to  provide  confidence  in  GOES  in-orbit  calibration  procedures. 

The  model  was  used  to  examine  the  responsivity  of  two  SN03  instrument  detectors. 
Results  indicate  that  the  percent  non-linearity  of  the  response  of  these  detectors  and  their  pre¬ 
amplifier  circuits  can  reach  4-6%,  which  is  significant  to  calibration.  Both  the  detector  and  the 
pre-amplifier  circuit  contribute,  and  contribute  additively,  to  the  percent  non-linearity.  The 
detector  contribution  was  dominant,  being  as  much  as  =70%. 

It  was  found  that  an  OGE  calibration  equation  of  quadratic  order  suitably  describes  the 
non-linearity  due  to  HgCdTe  detectors  and  detector  preamplifer  circuitry.  The  calibration 
equation's  quadratic  coefficient  was  found,  itself,  to  have  up  to  a  quadratic  dependence  on 
background  power  and  detector  temperature,  indicating  that  a  quadratic  interpolation  technique 
would  be  suitable  for  its  in-orbit  adjustment  to  variation  in  these  parameters. 
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APPENDIX  A 


HgCdTe  DETECTOR  MODEL 

This  appendix  summarizes  the  detector  physics  and  computational  procedures  used  for 
calculating  GOES  detector  and  detector-circuit  responsivity.  It  duplicates  some  of  the  discussion 
in  the  main  body  of  the  report.  It  is  intended  as  a  bridge  between  die  main  body  of  the  report  and 
the  actual  code  employed  given  in  Appendix  B. 

1  ENERGY  GAP,  QUANTUM  EFFICIENCY,  AND  Cd  MOLAR 

FRACTION 

The  bandgap  energy,  Eg,  vs  the  temperature,  T,  in  degrees  K  and  molar  fraction  of  Cd,  x, 
is  taken  from  the  work  of  Hansen, Schmit,  and  Casselman[Ref  1]. 

%  =  (-0.302  +  0.000535  T  (1  -  2  x)  +  1.93  x  -  0.81  x2  +  0.832  x3 ),  eV  (Al) 


The  gap  energy  equals  the  energy  of  a  photon  at  a  "cut-off'  wavelength,  Xc0  given  by 

^0=1*2222., 

h 


(A2) 


The  quantum  efficiency,  T),  gives  the  ratio  of  photo-excited  carriers  to  the  photon  flux 
intensity  impinging  on  the  detector.  The  expression  used  is  taken  from  the  report  by  Schoolar, 
BickneU,  and  Spears  and  references  cited  therein[Ref  2]. 

■n  =  (1  R)  l~Exp[-2aH]  .  ^ 

^  1  }  1  -  RExp[-2aH]  (  5) 


where 


a  =  absorption  coefficient,  cm*l 
H  =  detector  height  =  0.0008  cm 
R  =  front  surface  reflection  coefficient 

The  value  for  the  absorption  coefficient  at  a  wavelength  for  which  the  photon  energy 
equals  Ep  is  taken  from  work  reported  by  Schoolar  in  Reference  2. 


cto(l+15V^)-Eg  J,  cm-1 

(A4a) 

Exf^] 

,  cm-1 

Ep<Eg 

(A4b) 
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where: 


Oq  =  177(WEg,  cm*  1 
Tin  degree K 
Eg  and  Ep  in  eV 

k  =  Boltzman’s  constant  =  1.38066  x  10*^3  Joule/degree  K 
1  Joule  =  1/(1.60218  x  10‘19)  eV 
k  =  8.617  x  10’5  eV/degree  K. 

Quantum  efficiency  determines  the  wavelength  dependence  of  the  responsivity  of  the 
detector.  Using  relative  spectral  responsivity  data  collected  by  Ill'  and  the  above  equations,  an 
estimate  of  x  has  been  obtained  for  two  GOES-NEXT  SN03  detectors.  The  estimate  was  made 
by  graphically  fitting  relative  responsivities  predicted  by  the  equations  at  two  operating 
temperatures  to  the  measured  50%  relative  responsivity  points  in  the  cut-off  wavelength  region. 
The  value  of  R  =  0  was  used  for  the  estimated  match.  The  estimates  are  listed  in  Table  Al. 
Values  of  x  used  in  Reference  2  in  describing  expected  GOES-NEXT  detector  performance  are 
also  listed. 

TABLE  Al 


Estimated  values  of  x 


Detector 

Temperatures(K) 

Estimated  x 

[Ref  2] 

MW-L002 

102, 80.9 

0.2120 

— 

LW-J033 

102, 84.3 

0.1956 

0.1913 

L006 

105, 83 

0.2338 

— 

I4-L002 

105, 84 

0.2080 

0.209 

I5-J019 

105, 82 

0.2035 

0.203 

The  50%  relative  response  points  agree  at  the  two  temperatures  to  about  ±  2%  in 
wavelength.  Although  indicated  to  4  decimal  places,  the  accuracy  of  the  estimated  value  of  x  is 
no  better  than  ±  5%. 

2  INTRINSIC  AND  EQUILIBRIUM  CARRIER  CONCENTRATION 

Intrinsic  carrier  concentration,  nj,  is  taken  from  the  work  of  Hansen  and  Schmit[Ref  3]. 


ai  =  5.585  +  0.001753  T  -  3.82(x)  -  0.001364  (xT) 


bi  =  1014 


EgO-75  TlS 
Exp[5802.84  5? 


nj  =  ajbj, 


# 

cm3 


(A5) 


Equilibrium  density  of  electrons,  n0,  and  holes,  p0,  in  the  conduction  and  valence  bands 
in  the  case  of  a  material  donor  impurity  density,  Nd,  is  given  by 
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The  effective  carrier  relaxation  time,!,  is  determined  by  both  Auger  recombination  and 
surface  recombination  processes. 

4.1  AUGER  RECOMBINATION 

The  dominant  mechanism  in  carrier  recombination  is  the  Auger  process[Ref  5].  The 
intrinsic  population(no=Po=ni)  Auger  relaxation  time,  t^i»  *s  given  by 

,  Sec  (A9a) 

where  the  dependence  of  x^i on  the  temperature,  bandgap  energy,  and  electron  effective  mass  is 
taken  from  Reference  5.  The  effective  mass  of  a  conduction  band  electron  is  taken  to  be 
proportional  to  the  bandgap  energy  and  a  scaling  coefficient  has  been  chosen  to  give  agreement 
with  data  in  Reference  5  using  the  above  Eqs.(Al)  and  (A5)  for  nj[x,T],  and  Eg[x,T][Ref  6]. 

For  the  non-intrinsic  case  where  np  is  the  electron  carrier  concentration  induced  by 
photon  flux  the  Auger  relaxation  time,  xa,  is  given  by 
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tA-7 - ^"i^Ai  - 2WAl  ,  Sec  (A9b) 

(no+np](np+Po+np)  (no+npXno+Po+np) 

4.2  SURFACE  RECOMBINATION 

As  described  in  Reference  2  and  references  cited  therein,  surface  recombination  alters 
carrier  recombination  lifetime.  The  effect  of  surface  recombination  depends  on  the  ambipolar 
diffusion  constant,  Da,  given  by[Ref  7] 


Da  =8.617  10~5t 

(no+np)Pe+{Po+np^lh 


cm 

’Sec 


(A10) 


where 


^  =  8.617  10*5  T,  foyk— ,  T  in  degrees  K. 
q  Coulomb 

A  parameter  that  depends  on  Da  is  the  ambipolar  drift  length,  La,  given  for  a  earner 
recombination  time,  ta,  by 

La  =  VD^X,  cm  (All) 

The  effect  of  surface  recombination  on  the  effective  carrier  relaxation  time,t,  is  given  as 

t  =  Ajta  (A12) 

where 


Sinh 

H 

].StA 

La 

Cosh 

IL' 

=La= 

-1 

H 

(StA 

lLa 

n 

Sinh 

+  2Sta  Cos! 
La 

5  FLUX  DEPENDENT  CARRIER  LIFETIME 

The  flux  dependence  of  the  carrier  lifetime,  x,  is  the  principal  cause  of  photo-conductor 
response  non-linearity. 

5.1  FLUX  DEPENDENT  AUGER  LIFETIME 

As  indicated  by  Eq.(A9b)  and  as  experimentally  determined  by  Bartoli  et  al,  the  Auger 
lifetime,  ta>  and  thus  the  effective  carrier  lifetime,  t,  depends  on  the  carrier  concentration 
generated  by  photon  flux,  np[Ref  8].  For  an  incident  photon  flux  Qp,  #/cm2-Sec,  the  value  of  np 
is  given  by 
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(A14) 


=  ]ltQp 

p  H  ’  cm3 

Combining  Eq.(A9b)  with  Eq.(A12)  and  substituting  the  result  into  Equation  (A14)  leads  to  the 
following  cubic  equation  for  np: 


where 


tAi=/cAiAi  • 

The  solution  of  Eq.(A15)  satisfies  both  Equation  (A9b)  and  (A  14).  The  effect  of  surface 
recombination  introduces  two  computational  complications:  First,  equation  (A  15)  must  be 

solved  to  evaluate  ta  from  Eq.(A9b).  But,  ta  must  be  known  in  order  to  evaluate  La  from 
Eq.(All)  and  A\  from  Eq.(A13).  A  second  complication  is  that  La  depends  on  Da  which,  as 
indicated  in  Eq.(A10),  depends  on  np. 

A  solution  to  this  problem  was  accomplished  by  Marion  Reine  of  (now)  LIRIS  in  1979 
for  ATMOS  detectors  operated  at  77K[Ref  9].  Marion  used  a  different  approach  for  accounting 
for  the  surface  recombination  effect,  but  incorporated  its  flux  dependence. 

S2  ITERATIVE  SOLUTION  FOR  x  AND  np 

The  following  is  the  iterative  technique  used  for  self-consistently  solving  Eqs.  (A9b), 
(A10),  (All),  (13),  and  (A15): 

Let  np(j]  be  the  jth  estimate  of  np. 

(i)  Given  that  {no.po^Ail  316  calculated  for  a  given  {Nd,x,T},  an  estimate  of  ta,  taUL 
may  be  obtained  from  Eq.(A9b)  using  np[j]. 

(ii)  Also  using  np[j],Eq.(A10)  may  be  evaluated  to  give  an  estimate  of  Da,Da[j].  Then, 
using  Da[j]  and  taOL  Eq(Al  1)  may  be  solved  for  an  estimate  of  La,  La[j]. 

(iii)  Using  La[j]  and  TaU]>  Eq(A13)  may  be  solved  for  an  estimate  of  Aj,  Ajjj]. 

(iv)  The  value  of  Ajjj]  may  then  be  used  to  compute 

=  tAi  A{j]. 


(v)  Letting  Np  { n0,po>Qp,t* AiUl  1  denote  the  single  positive  np  root  of  Eq.(A15), 
the  (j+l)th  estimate  of  np  is  given  by 

np[j+l]  =  Np{no,po,Qp/t*AiIj]} 
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and  the  iteration  is  repeated. 


The  goal  of  the  iteration  is  to  accurately  estimate  z  for  a  given  Qp  and  to  accordingly 
assess  np.  Upon  the  iterative  procedure  converging,  the  iteration  is  ended  when 


tQ+1]  i 
T[j] 


=  1%. 


The  initial  value  of  np,  np[0]  is  taken  from  Eq  (A14)  with  x[0]  =  t^i-  A  solution  is 
described  by  the  order  of  the  iteration  undertaken.  Results  are  reported  using  10th  or  higher 
order  solutions. 


6  DETECTOR  REPONSIVITY 

Determining  the  responsivity  of  the  GOES  detectors  requires  separating  Qp  into  two 
components, 

Qp  =  Qb  +  Qs 

where  Q5  is  a  background  flux  that  is  assumed  constant  on  the  detector  both  between  and 
throughout  measurements  and  Qs  is  the  signal  flux  incident  on  the  detector  during  a  look  at  an 
earth  scene.  Thus,  from  Eq.(A14)  the  photo-excited  carrier  concentration,  np,  is  also  composed 
of  two  components,  n^  and  ns  from  the  background  and  the  signal  flux,  respectively. 

The  relation  between  incident  flux,  Qp,  and  the  incident  power,  Pp,  in  watts  at  the  cutoff 
wavelength,  Aco>  is 


=  Wp  # 

VP  hcWL’  cm2-Sec 


(A16) 


where 

h  =  Planck’s  Constant  =  6.63x  10' 34  Joule-Sec 
c  =  speed  of  light  =  3x  10*4 

and  where  W  is  the  width  and  L  the  optical  sensitive  length  of  the  detector.  Thus,  the  power 
incident  on  the  detector  is  composed  of  two  components,  Pb  and  Ps,  from  the  background  and 
the  signal,  respectively. 

The  responsivity  is  evaluated  at  the  peak  response  wavelength  which  is  assumed  here  to 
be  Xqq.  The  responsivity  is  evaluated  in  two  steps: 

First,  a  voltage  is  computed  for  np  =  nj).  This  voltage  will  depend  on  the  detector  resistance, 
Rd[np]  for  an  incident  flux.  R<j[npJ  is  given  by: 
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Rd["pl  J2NaHss 
ql  H 


_J _ . 

+Mno + np)l^- 


Ohms 


(A  17) 


where 


q  =  1.60218x10*19  Coulombs 
H  =  detector  width,  approximately  =  0.8x10*3  cm 
W  =  detector  width 

L  =  detector  photo-responsive  length  =W  =  (5.0-7.5)  10*3  cm 

W*L  =  photo-responsive  area 

Le  =  detector  electrical  length,  =(2-4)L 

(NSs  ms) =  Surface  State  (population-density  mobility)  product=(2-5)10l6,  #/Volt-Sec 

The  value  of  the  (Nss  |iss)  product  depends  on  propietary  fabrication  processes  and  the 
estimate  given  is  only  an  edjucated  guess.  The  estimated  values  give  a  shunt  conductance  value, 
gss,  where 


Sss  —  2  q(Nss  pss)W/Le, 


of  approximately  (4-7)  mMho  for  SN03  detectors.  Values  are  estimated  from  detector  resistance 
vs  bias  and  temperature  data. 


Given  R<j[np]  the  detector  preamplifier  voltage, V[Pp],  is  determined  by  the  bias  current, 
Iq,  the  bias  voltage,  Vc,  and  the  bridge  bias  resistor  value,  R^. 


V[Pp]  =  I0 


RdIPp] 
Rb  +  RdP  pi 


-,  Volts 


(A  18) 


Io  =  ^, Amps  and Vc  =  10 Volts. 

Rb 

The  second  step  is  the  evaluation  of  V[Pb  +  Psl,  which  follows  from  Eq.(18)  with  Pp  =  (P5  + ' 
Ps).  The  preamplifier  earth  scene  signal  voltage,  AV[Psl,  is  then  evaluated  as  the  difference 
between  V[Pb  +  Psl  and  V[Pbl- 


AV[PS]  =  V[Pb  +  Ps]  -  V[Pb].(19)  (A19) 

Equation  (A  19)  gives  the  basic  expression  for  the  investigated  GOES  SN03  responsivity 
non-linearity. 

7  DETECTOR  SMALL-SIGNAL  RESPONSIVITY  AND  PC  GAIN 

Two  parameters  that  are  of  interest  in  assessing  detector  performance  are  the  small-signal 
responsivity  and  the  photo-conductive(PC)  gain. 

A  small-signal  voltage  responsivity,  5R[Pb,0],  is  taken  to  be 
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l  imit  AV  dV[Pp]  dV[npl  (tap  dQp  _  T)XC0  dVtnp]  d(Qpt[Qp| 

9i[Pb,0]-  p  Ps  ”  dPp  "  (tap  dQpdPp  hcHWL  dnp  dQp  I  Qp^Qb 

P  (A20) 

Equation  (A20)  indicates  that  the  derivative  of  t[Qp]  evaluated  at  Qp  =  Qb  plays  a  role  in  the 
small-signal  responsivity,  which  is  bom  out  by  modeling  results.  Taking  Rb»R(j[Pb] 
Eq.(A18),  the  small-signal  responsivity  is  given  by  Eq.(A20)  to  be 

StPb.O]  =  -0.8  0  Io  Me  %ffPb]  rP  (MPbl)2.  ^  (A21) 

6 

where 

XefftQb]  =  -IQbl-Qb^i| 

In  the  computer  code,  the  small  signal  gain  for  a  given  Pb  is  evaluated  by  determining  the 
responsivity  for  a  signal  power  that  is  0.00000001%  of  the  background  power.  This  small  signal 
flux  responsivity  is  taken  to  be  the  responsivity  the  detector  approaches  as  Ps->  0. 

The  photo-conductive  gain,  GpC,  is  defined  as  the  ratio  of  the  mean  free  path  of  a  photo- 
excited  carrier  to  the  photo-electrically  sensitive  length(L)  of  the  detector.  When  a  carrier 
transits  the  detector  and  is  swept  into  a  collecting  contact,  charge  conservation  requires  that  the 
opposite  contact  emit  a  carrier.  The  idea  is  that  a  single  photo-electrically  excited  earner  may 
effectively  transit  the  detector  several  times.  The  gain  gives  the  number  of  transits  for  L=Le- 

The  mean  free  path  of  the  carrier  is  just  its  effective  lifetime,  teff,  multiplied  by  its  mean 

velocity,  peE,  where  E  is  the  electric  field  due  to  the  bias  voltage.  In  this  case  the  voltage  across 
the  detector  is  given  by  the  product  (IqR^)  leading  to 

0pc=i(^pVfiPbl) 

So,  considering  Eq.(A21) 

9t[Pb.O]  s  -0.8XCO  T1  GpC  R<i  , 

Measurement  of  the  responsivity  and  the  detector  resistance,  and  hence  the  inference  of  the  (T| 
Gpc)  product,  is  often  undertaken  to  characterize  detectors. 
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APPENDIX  B 


COMPUTER  CODE 

Calculations  employed  the  code  listed  in  this  Appendix.  The  code  is  written  for 
execution  in  Mathematica  and  is  similar,  but  not  identitical  to,  C  language.  Mathematica  can 
convert  this  code  to  C  or  to  Fortran. 

It  is  suggested  that  anyone  wishing  to  calculate  detector  performance  compose  their  own 
code.  This  listing  is  given  to  serve  as  a  guide  to  one  way  it  can  be  done. 

The  code  is  comprised  of  two  functions:  fO  and  perfP[].  The  functions  take  as  input  the 
variables,  listed  within  the  brackets!].  Variables  have  to  start  with  a  lower  case  character  and 
Greek  symbols  and  subscripts  are  not  supported  by  Mathematica  code.  Thus,  Nd->nD,  S->s, 

Qp->qp,  eg->eg,  X->wvl,  etc. 


The  first  function,  f[|,  takes  basic  the  parameters  [x,  T,  Ntf,  S,  Qp]  and  solves  for  the 
equilibrium  carrier  population  density,  np  and  the  carrier  lifetime.  It  does  so  using  a  Do-loop 
that  cycles  ”kay"  times.  In  practice  I  have  found  kay=10  gives  satisfactory  results.  The  function 
returns  these  values  and  other  intermediate  results  that  are  useful  in  calculating  the  responsivity 
of  the  detector. 

f[x_,t_,nD_,s_,qp_,kay J  := 

( *Define  variables  used  within  the  function*) 

Module!  { eg,wvl,ai,bi,ni,no,po,tAi,muE,muH,a,b,np,tAug,do,doN,doD,la,xa,xb,aA,tAiE,c  joot } , 

(*Calculate  fixed  Quantities* ) 

eg=-0.302+1.93*x-0.810*(xA2)+0.832*(xA3)+5.35*(10A-4)*t*(l-2*x); 
wvl= 1.2399/eg; 

ai=5.585-3.82*x+(l  .753*  10A-3)*t-(l  .364* 10A-3)*x*t; 

bi=(10A14)*(egA0.75)*(tA1.5)*Exp[-eg/(2*(8.617*10A-5)*t)]; 

ni=ai*bi; 

no=(nD/2)+Sqrt[((nD/2)A2)+niA2]; 
po=-(nD/2)+Sqrt[((nD/2)  A2)+niA2] ; 

tAi=(3*10A-13)*(((eg/((8.617*10A-5)*t))A1.5)/eg)*Exp[eg/((8.617*10A-5)*t)]; 

muE=(2.5*(10A-7)*(xA2)*(tA1.5)+3.13*(10A-2)*(xA4)*(tA-1.5))A-l; 

muH=(5.26*(10A-7)*(xA-l)*(tA1.5)+4.0*(10A3)*(xA7)*(tA-1.5))A-l; 

(* All  fixed  quantities  are  calculated->Next  enter  into  iterative  calculation *) 

(*Do-loop  statement*) 

Do[ 

a=N[(2*no+po)/ni] ; 
b=N[  1 +(no/ni)A2] ; 

np[0]=  ((2*tAi*(niA2))/(no*(no+po)))*(qp)/0.0008; 
aA[0]=l; 

tAug=2*tAi*aA[k-l]*(niA2)*(((no+po+np[k-l])*(no+np[k-l]))A-l); 


41 


doN=8.617*(10A-5)*t*(no+po+2*(np[k-l]))*muH*muE; 

doD=((no+np[k-l])*muE)+((po+np[k-l])*muH); 

do=doN/doD; 

la=(do*tAug)A0.5; 

xa=0.0008/la; 

xb=s*tAug/la; 

aA[k]=N[(Sinh[xa]+xb*(Cosh[xa]-l))*(((l+(xbA2))*Sinh[xa]+2*xb*Cosh[xa])A-l)]; 

tAiE=tAi*aA[k]; 
c=N[(2*tAiE/(0.0008*ni))*qp]; 
root=  {ToRules[NRoots[yA3+a*yA2+b*y-c=0,y]] } ; 
np[k]=ni*(y/.root[[3]]),{k,  1  ,kay }]; 

(* Finish  of  Do-loop  statement *) 

tAugO=2*tAi*aA[kay]*(niA2)*(((no+po+np[kay])*(no+np[kay]))A-l); 

(* Return  list  of  results  to  12  place  precision*) 

N[  { tAugO,np[kay]  ,muE,wvl,no,ni,tAugO*qp/0.0008 } ,  1 2] 

] 


The  function  perfPQ  takes  as  an  argument  the  parameters  [x,  T,  Nd,  S,  Io,  gSs(surface 
state  shunt  conductance  in  Mhos),  Adet(area  of  det  in  cm^),  r(ratio  of  ho/L),  Pb,  P si 
calculates  the  voltage  generated  by  the  detector  for  incident  flux  due  to  Pb,  PbQ+lO^lO),  and 
Pb+Ps-  From  these  results  the  small  and  large  signal  responsivites  are  calculated.  The  percent 
change  from  small  to  large  signal  responsivity  is  assessed.  Also,  the  percent  change  in  the 
effective  carrier  lifetime  for  small  and  large  signal  response  is  evaluated.  The  function  returns  a 
list  of  several  parameters  that  are  of  interest  to  the  detector  and  detector  preamplifier  operation. 

perfP[x_,t_,nD_,s_,io_,gss_,aDet_,r_,pb_,ps_,kayJ:= 

(*Deftne  variables  used  within  the  function*) 

Module[{bl,b2jB^b^Db^DbsjDbSS,vb,fbs,vbs,sig,resp,respNL,percentNL,percenttau,qbSS,vb 
S  S  ,respS  S  ,qb,qs,wvl,eg } , 


rB=(10.0)/io; 


(*  Calculate  fixed  quantities*) 


eg=-0.302+1.93*x-0.810*(xA2)+0.832*(xA3)+5.35*(10A-4)*t*(l-2*x); 
wvl= 1.2399/eg; 


qb=pb*  wvl/(aDet*2*(  1  1 9)); 

qs=ps*wvl/(aDet*2*(10A-19)); 

(*CallfI]for  Qp<->Pb,  compute  Vb*) 

fb=f[x,t,nD,s,qbJcay]; 

rDb=(go+(l/r)*(1.6*10A-19)*(fb[[3]])*(fb[[5]]+fb[[2]])*0.0008)A-l; 

vb=10.0*rDb*((rB+rDb)A-l); 

(*Callf[]  for  Qp<->Pb+Ps,  compute  Vbs*) 

fbs=f[x,t,nD,s,qb+qs,kay] ; 

rDbs=(go+(l/r)*(1.6*10A-19)*(fbs[[3]])*(fbs[[5]]+fbs[[2]])*0.0008)A-l; 
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vbs=10.0*rDbs*((rB+rDbs)M); 


(* Evaluate  resp[Pb+Ps]*) 

sig=(vbs-vb)*  10*6; 

resp=(fb[[4]]*(vbs-vb))/(2.0*(10A-19)*aDet*qs); 

(*Callf[]  for  Qp<->(1 +10*-10)Pb, compute  VbSS*) 
qbSS=qb*(1.0000000001); 
fbSS=f[x,UiD,s,qbSS,kay]; 

rDbSS=(go+(l/r)*(1.6*10*-19)*(fbSS[[3]])*(fbSS[[5]]+fbSS[[2]])*0.0008)*-l; 

vbSS=10.0*rDbSS*((rB+rDbSS)A-l); 

(*Evaluate  small-signal  resp*) 
respSS=(fb[[4]]*(vbSS-vb))/(2.0*(10A-19)*aDet*((qbSS-qb))); 

(*Compute  quantities  of  interest*) 

bl=l/resp; 

percentNL==((resp/respSS)- 1  )*  1 00; 

b2=((bl)A3)*100*(2.0*(10M9)*aDet*qs)/(fb[[4]]*percentNL); 

percenttau=((£bs[[l]]/fb[[l]])-l)*100; 

(* Return  list  of  results  of  interest  to  6  place  precision*) 

N[  { sig,resp,respSS,rDb,rDbs,qb,qs,t,pb,ps,percentN,percenttau,bl,b2},6] 

] 


To  operate  the  code  it  is  necessary  that  both  f[]  and  perfP[]  are  evaluated,  i.e.  "loaded", 
into  the  Notebook  or  "routine"  being  run  to  do  the  calculations.  The  function  perfP[]  calls  for 
output  from  f[]. 
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APPENDIX  C 


ADDENDUM 

During  the  completion  of  this  report  data  from  preliminary  vacuum  tests  of  the  serial 
number  3(SN03)  Imager-3  detector  became  available.  The  data  were  collected  with  the  detector 
in  the  instrument.  The  instrument  responsivity  was  measured  vs  detector  temperature. 

This  was  of  particular  interest  since  the  model  predicted  that  the  Imager  3  responsivity 
could  increase  over  the  detector  operating  temperature  range  with  increasing  detector 
temperature.  (See  Figure  10.)  This  behavior  is  somewhat  counter-intuitive  to  the  notion  that 
HgCdTe  detector's  performance  improves  with  decreasing  temperature:  In  the  case  where 
detector  responsivity  dominates  performance,  the  data  indicate  that  this  is  not  necessarily  true,  at 
least  for  all  HgCdTe  detectors. 

A  plot  of  the  data(points)  and  the  model  fit(solid  line)  is  shown  in  Figure  C-l. 


Figure  C-l.  Fit  of  Imager  3  model  to  in-situ  instrument  test  data  of  normalized 
responsivity  vs  detector  temperature.  ( Assumed  background  power,  Pb  =21nW.) 

To  fit  the  model  the  value  of  the  impurity  density,  Nd(Table  1)  was  varied  from  the 
assumed  5xl0^/cm3  to  6.5xl0^/cm3.  It  can  be  seen  from  Figure  C-l  that  a  value  of 
=6xl014/cm3  gives  qualitative  agreement.  And  it  can  be  seen  that  the  responsivity  of  Imager-3 
does  increase  with  increasing  detector  temperature  over  this  range,  as  predicted  by  the  model. 

An  explanation  of  this  effect  is  given  on  pages  21-23  of  the  text.  A  consequence  of  the 
effect  is  that  lowering  in-orbit  detector  temperature  may  improve  Sounder  instrument 
performance  in  the  Long-Wavelength  bands,  but  may  also  slightly  degrade  Imager-3 
performance. 

The  author  appreciates  Richard  Heinrichs’  providing  a  file  of  the  test  data. 
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